Draft version January 6, 2012 

Preprint typeset using IAT^X style emulateapj v. 5/25/10 



ON THE NATURE OF THE TRANSITION DISK AROUND LKCA 15 

Andrea Isella, Laura M. Perez, John M. Carpenter 

Department of Astronomy, California Institute of Technology, MC 249-17, Pasadena, CA 91125, USA 

Draft version January 6, 2012 

ABSTRACT 

We present CARMA 1.3 mm continuum observations of the T Tauri star LkCa 15, which resolve 
the circumstellar dust continuum emission on angular scales between 0.2 ,, -3 ,/ , corresponding to 28-420 
AU at the distance of the star. The observations resolve the inner gap in the dust emission and reveal 
an asymmetric dust distribution in the outer disk. By comparing the observations with theoretical 
disk models, we calculate that 90% of the dust emission arises from an azimuthally symmetric ring 
that contains about 5 x 10 -4 M of dust. A low surface-brightness tail that extends to the northwest 
out to a radius of about 300 AU contains the remaining 10% of the observed continuum emission. 
The ring is modeled with a rather flat surface density profile between 40 and 120 AU, while the inner 
cavity is consistent with either a sharp drop of the 1.3 mm dust optical depth at about 42 AU or a 
smooth inward decrease between 3 and 85 AU. We show that early science ALMA observations will 
be able to disentangle these two scenarios. Within 40 AU, the observations constrain the amount of 
dust between 10 -6 and 7 Earth masses, where the minimum and maximum limits are set by the near- 
infrared SED modeling and by the millimeter-wave observations of the dust emission respectively. In 
addition, we confirm the discrepancy in the outer disk radius inferred from the dust and gas, which 
corresponds to 150 AU and 900 AU respectively. We cannot reconcile this difference by adopting an 
exponentially tapered surface density profile as suggested for other systems, but we instead suggest 
that the gas surface density in the outer disk decreases less steeply than that predicted by model fits 
to the dust continuum emission. The lack of continuum emission at radii lager than 120 AU suggests 
a drop of at least a factor of 5 in the dust-to-gas ratio, or in the dust opacity. We show that a sharp 
dust opacity drop of this magnitude is consistent with a radial variation of the grain size distribution 
as predicted by existing grain growth models. 



1. INTRODUCTION 

Over the past decade, the Spitzer Space Telescope 
has identified an increasing population of circumstellar 
disks surrounding pre-main sequence stars that display 
a deficit of near- and mid-infrared flux compared to a 
"classical" T Tauri star, while exhibiting similar levels 
of far-infrared emission. The spectral energy distribu- 
tion (SED) of these objects suggests that these stars lack 
warm dust in the inner disks regions, while substantial 
emission at longer wavelengths indica tes that the oute r 
disk remains quite massive (see, e.g., Calvet et al.|2 005). 
These so-called transition disks may be intermediate sys- 
tems between young (~ 1 — 10 Myr), optically thick, gas 
rich protoplanetary disks and old (> 10 Myr), optically 
thin, gas poor, debris disks, where significant disk evo- 
lution has occurred and a planetary system might have 
formed. 

Subsequent high-resolution images of the optically thin 
millimeter-wave dust emission confirm the SED model- 
ing predictions tha t the inner disk in transition disks 
is depleted of dust ft Andrews et al. 2011a ; Brown et al. 
2009; |Hughes et ai.j]2UU9| |Pietu et al.||2006| ). However, 
spatially resolved observations ol the gas emission show 
that some of these d ust depleted cavitie s contain a signif- 
icant amount of gas ( |Pietu et~ al. 2007) and that the size 
of the hole is not alwa ys consistent with SED modeling 
( Andrews et al.pOlla ). Furthermore, recent high angu- 
lar resolution observations of the mm- wave dust emission 
have also revealed that some "classical" disks have dust 



depleted inner regions similar to those of "transition" 
disks , but do not show any deficit in the infr ared ex- 
cess fllsella et aL]|2010a|b| [Andrews et al.||2011a| ). These 
results demonstrate that SEDs alone cannot provide a 
complete picture of the disk structure. High angular res- 
olution millimeter-wave observations are essential to re- 
move some of the degeneracies that limit SED modeling, 
and therefore provide a more accurate description of the 
radial distribution of circumstellar material. 

To determine the physical phenomena responsible for 
disk cavities, measurements of the mass surface den- 
sity profile E(r) are essential. To first order, circum- 
stellar disks evolve via viscous spreading, where friction 
forces redistribute t he material and smoothly sha pe the 
disk surface density ( |Lynden-Bell fc Pringle||l974| ) . Disk 
photoevaporation by energetic photons might reduce the 
amount of disk material. Models that combine these two 
process have been the focus of many theoretical studies 
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(I Alexander et al.||2006a|b| |Gorti et al.|2009||Owen et al. 
|z010[ ). These models suggest that viscous evolution pro- 
ceeds unperturbed by photoevaporation until the mass 
accretion rate throughout the disk is comparable with 
the photoevaporation rate. After this phase, which can 
last for a few Myr for typical values of the stellar radia- 
tion and disk surface density, a gap opens at few AU from 
the central star, where the photoevaporation rate is the 
largest. The material within the gap then drains quickly 
onto the central star leaving a gas-free cavity in the disk. 
Although the exact time scale for the formation of the 
cavity varies between different models, they all predict 
an order of magnitude drop in the gas surface density 
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over a 1 AU radial extent at the edge of the cavity. 

Alternatively, cavities in the disk can be produced by 
the gravitational interacti on between giant p lanets and 
the disk material (see, e.g., Bryden et al.|1999| ). Alth ough 
this i s a co mmon explanation lor transition disks, |Zhu| 
et al.| ( |2Q11[ ) have shown that even multiple giant planets 
cannot account for optically thin cavities larger than 20 
AU, unless the small grains within the region perturbed 
by the planets are efficiently removed. Similarly to the 
photoevaporation model, the planet-disk interaction pre- 
dict a sharp drop of the dusty and gas density at the edge 
of the region cleared by the planet tidal interaction. 

A third possible origin of the inner disk cavities is that 
the dust grains have grown from micron sized particles to 
centimeter sized pebbles. Therefore, the cavities reflect a 
reduced dust opacity and not a reduced matter density. 
Theory predicts that dus t coagulation is more efficient 
in the dense inner disk ( jTanaka et al 



2005 



Dullemond 



Dominik||2005 ) and as the dust particles grow, the 
opacity ol the dust diminishes. Current models of grain 
growth anticipate a continuous r adial variation of th e 
dust opacity within the disk (e.g., [Birnstiel et aT]|2010[ ), 
whose observational signature corresponds to a gradual 
transition between the inner and outer disk structure. 

Among the known stars with a transition disk, the 2- 
5 Myr old star LkCa 15 represents an outstanding case. 
LkCa 15 is a K5 star (L* = 0.741^, = l.OM ; |Sirnon| 



et al.||2000 
Taurus star- 



Kenyon fc Hart mannj [1995]) located in the 



orming region at a distance of about 140 pc 



van den Ancker et al.|1998[ ). Its disk is characterized by 



a parti ally dust-d epleted inn er region of about 50 AU in 
radius (|Andrews et al.|2011a Espaillat et al.|2007 ; Pietu 



et al.|[2QQ6[) and a mass accretion rate of about 3 x 10 y 
M yr" 1 (Hartmann et al.|l998[ ). The disk SED is char- 
acterized by the presence ol an infrared excess over the 



y tne presence 
stellar photosphere, which can be explained by the pres- 
ence of hot dust with i n few AUs from the central star 
(lEspaillat et all 120071 iMulders et al.||2010|). The disk 



has a very acti ve chemistry observed m several molecu- 



lar transitions ([Chapillon et al.||2008 Pietu et al. 2007 



Oberg et al. 2010), and a lso exibits keple rian rotation 



out to a radius of 900 AU flPietu et al.l|2007| ). These lat 



ter observations suggest that gas is present at an orbital 
radius as small as 10 AU, and therefore well within the 
cavity observed in the dust continuum emission. Ground 
based H-band (1.6 fim) observations reveal the presence 
of an elliptical structure i n correspondence to the edge 
of the dust inner cavity (Thalmann et al. 2010). Fi- 



nally, Kraus fc Ireland| ( |201ip has recently reported the 
discovery through aperture masking of a possible 6 Mj 
mass proto-planet located inside the continuum cavity 
at an orbital distance of 16 AU from the central star. 
These observations suggest that the cavity might be in- 
deed shaped by the dynamical in tera ction with a giant 
planet. Howev er, |Zhu et aT| ( |20 11 [ ) and |Dodson-Robinson| 
Salyk ( |2011[ ) have argued that a cavity 50 A U in radius 
cannot be explained by a single planet, and that addi- 
tional giant planets or clearing mechanisms are required 
to explain these observations. 

In this paper we present new observations of the 1.3 
mm dust thermal emission obtained with the Com- 
bined Array for Research in Millimeter-wave Astronomy 
(CARMA). The observations achieve an angular resolu- 



tion of 0.2" corresponding to a spatial scale of 28 AU 
at the distance of star, which is a factor of two improve- 
ment with respect of existing observation of this system. 
In Section[2]and|3]we present our observations of LkCa 15 
and discuss the morphology of the dust emission. The 
method used to analyze the data is discussed in Section[4j 
and the results are presented in[5j In Section|6]we discuss 
the possible origin of the dust continuum cavity and of 
the whole disk structure. A summary of the main results 
follows in Section 

2. OBSERVATIONS AND DATA REDUCTION 

LkCa 15 was observed with CARMA array configu- 
rations A, B, C, and E, providing projected baseline 
lengths between 6 and 1700 m. C-configuration observa- 
tions were obtaine d in October 2007 and are discussed in 
|Isella et aL| ([2009 ) . Observations in B-conflguration were 
obtained in December 2008. The receivers were tuned 
at a local oscillator frequency (LO) of 224.75 GHz and 
the continuum emission was observed using three con- 
tiguous bands of 468 MHz each in both correlator side 
bands, resulting in a total bandwidth of about 3 GHz. 
Flux and pass band calibration were derived by observing 
Uranus and 3C84 respectively; 0510+180 was observed 
every 15 minutes to derive the atmospheric and instru- 
mental phase corrections. A- and E-configuration obser- 
vations were carried out in November 2010 and January 
2011 respectively, using the upgraded CARMA correla- 
tor which provides a maximum continuum bandwidth of 
about 8 GHz. The receivers were tuned to the same LO 
frequency used for C and B-configuration observations. 
The upper and lower correlator side bands were config- 
ured with eight partially overlapping bands of 494 MHz 
each in order to cover the full bandwidth. The band pass 
shape was calibrated by observing 3C84 and 3C454.3; 
flux calibration was set by observing Uranus and cross 
checked using almost simultaneous observations of 3C84 
and 3C454.3 obtained with the SMA at the frequency 
of 225 GHz. The estimated uncertainty in the absolute 
flux calibration is 10%. Atmospheric and instrumental 
effects were corrected by observing the nearby quasar 
0510+180 every 6-10 minutes, depending on the array 
configuration. 

The data taken in the four different array configura- 
tions were independently calibrated using the MIRIAD 
software package. The dataset were shifted both in right 
ascension and declination to compensate for the stellar 
proper motion (see Appendix) and combined to increase 
the image sensitivity and dynamic range. Images were 
then obtained through standard inversion of the uv data 
using different weighting functions to emphasize the mor- 
phology of the dust emission on both large and small 
angular scales. The images were deconvolved using the 
CLEAN algorithm. 

3. MORPHOLOGY OF THE CONTINUUM EMISSION 

The high dynamic range obtained by combining the 
compact and extended array configurations enables 
CARMA to reveal the structure of the dust emission on 
both small and large angular scales. Figure [I] shows the 
map obtained by adopting three different weighting func- 
tions to achieve an angular resolution of (a) 3.0" x 2.6" 
(b) l.l"x0.5" and (c) 0.21"x0.19". In panel (a), the disk 
emission appears fairly symmetric and centered near the 
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Figure 1. Maps of the 1.3 mm dust continuum emission observed toward LkCa 15 employing three different weighting schemes of the uv 
data, (a) Natural weighting was adopted to achieve a FWHM beam size of 3.0 // x2.6 // . Intensity contours start at 3a and are spaced by 6<r, 
where the la noise level is 0.58 mJy/beam. (b) Robust weighting was adopted to achieve a resolution of l.l^xO.5". Contours are plotted 
every 3cr, where the lcr noise level if 0.53 mJy/beam. (c) Uniform weighting was adopted to achieve a FWHM beam size of 0.21" x0.19 /x . 
Contours are plotted every 3cr, where the la noise level is 0.43 mJy/beam. (d) The large, medium, and small red circles show the position 
of the center of symmetry of the dust emission corresponding to panel (a), (b), and (c) respectively. The symbol ★ marks the position of 
the star and the surrounding circle its uncertainty, calculated as discussed in the Appendix. The dotted curves show the intensity contours 
as in panel (c). 

position of the star. The integrated flux measured in 
a circular aperture of 4" in radius, is 128 ±5 mJy, and 



i t is consistent with the value measured by Pietu et al. 
(2006). In panel (b), the emission has two peaks sepa- 
rated by about 0.3". In addition, the low intensity con- 
tours show that the disk extends toward the northwest. 
The integrated flux, as measured in a circular aperture 
of 2" in radius is 127±15 mJy, where the uncertainty ac- 
counts only for random noise. Finally, the map shown in 
panel (c) reveals that the dust emission comes from an 
elliptical ring that extends for about 1.8" along its major 
axis, and for about 1.2" in the perpendicular direction. 
Assuming that the emission arises from a circular disk, 
the aspect ratio of the contour levels leads to a disk in- 
clination of about 50°, a position angle measured east 
from north of about 60°, and a disk diameter of about 
250 AU for the stellar distance of 140 pc. The total 
integrated flux measured in an elliptical aperture with 
major axis of 2.4" and the same aspect ratio of the disk 



is 117±10 mJy. CARMA observations reveal for the first 
time an asymmetric dust distribution in LkCa 15 outer 
disk, while the ring-like structure presented in panel (c) 
of Figure [T] is consistent with earlier lower angular resolu- 
tion observations of LkCal5 (Pietu et al. 2006; Andrews] 
et al.||2011aj ). 
lb mvestigc 



gate the properties of the dust emission we 
define the center of symmetry (x c ,y c ) of the map inten- 
sity I(x,y) such that x c = Y,xI/Y,I and y c = Ei/J/XJ, 
where the sum extends over all positions brighter than 5 
times the noise level in the image. Since I(x, y) depends 
on the weighting scheme adopted in the image formation 
process, the position of the center of the emission is a 
function of the angular resolution in the image. Panel 
(d) in Figure [I] shows the center of symmetry for the 
three maps discussed above. The stellar position and its 
uncertainty, as expected from proper motion correction 
(Appendix), are shown by the star symbol and the sur- 
rounding circle respectively. We find that the center of 



Isella A., Perez L., and Carpenter J. 



symmetry of the highest angular resolution map is consis- 
tent with the position of the star, and both appear to be 
shifted by about 0.1" with respect to the center of inner- 
most intensity contours levels. Although the astrometric 
uncertainty on the position of the star does not allow 
us to firmly conclude on the nature of this shift, we note 
that a pericenter offset of similar magnitude (0.064") and 
direction has been recently derive d from infrared obser- 
vations by Thalmann et al. (2010). 

At lower angular resolution the center of symmetry 
moves progressively toward the northwest up to a dis- 
tance of about 0.3" from the star. This shift of the cen- 
ter of symmetry, which results from the asymmetry ob- 
served at intermediate angular resolution, suggests that 
the outermost disk regions might have been perturbed. 
However, about 90% of the dust emission arises from the 
fairly symmetric ring shown in panel (c). The structure 
of this ring and the significance of the observed asymme- 
tries are discussed in the following two sections. 

4. DISK MODEL 

The disk properties were inferred using two prescrip- 
tions for the disk surface density. 

The first is a classical power law £(r) = S (ro/r) p 
characterized by four free parameters: the inner radius 
ri n , the outer radius r ou t: the normalization value Ho, 
and the slope p. This parameterization provides the 
simplest form for the disk surface density that mimics 
the sharp inner disk truncation caused by disk photo- 
evaporation or by dynamical interaction with a giant 
planet, but has the disadvantage of introducing an un- 
physical sharp edge of £(r) at the disk outer radius. 

To obtain a smoother form of S(r) we adopt a sec- 
ond parameterization given by the similarity solution 
for the disk surface density of a viscous keplerian dis k 
(see |Lynden-Bell fc Pringle|197 4; Hartmann et al. 1998), 
adopting the form presented in Isella et al. J 2009D, 



S(r) 



x exp • 



2(2-7) 



n 



(2-7) 



(!) 

This functional form has five free parameters, i.e., the 
disk inner and outer radius (r$ n and r ou t), the character- 
istic radius r tl the normalization value Yt t — E(r^), and 
the slope 7, which is related to the radial viscosity profile 
through the equation v(r) oc r -7 . 

Since the limited angular resolution and sensitivity of 
existing observations does not generally allow to fit for 
all five parameters, previous studies of the dust distribu- 
tion in transitional disks based on this form of th e sur- 
face density adopted two different approaches. In [isella 
|eTaD([200m|2010a|b|), we fixed the disk inner radius at 
the dust evaporation radius (< 0.5 AU), and compared 
the model to the observations to constrain the values of 
7, rt, and T, t . In this way, we found that the cavities 
in the dust continuum emission observed in LkCa 15, 
RY Tau, and MWC 758 disks might be explained with 
small inner disk radii and negative values of 7, which 
leads to a smoothly increasing surface density up to a 
radius r max = r + x (— 2 r y) 1 ^ 2 ~ ry \ A d ifferent approach 
was adopted by |Andrews et ah] ( |2011a| ), who fixed the 
slope 7 = 1, and derived the disk inner radius n n , r t and 



and The first modeling procedure, i.e., variable 7 
and fixed r^ n , stems from the idea that transitional disks 
might be shaped a global physical process, such as the 
disk viscosity and/or the variation of the dust opacity 
due to the grain growth. By contrast, the second case 
assumes that the disk surface density p rofile follows the 



a- const ant prescription of the viscosity ([Shakura 
yaev 19731) which leads to 7 



Sun- 

and that the cavities 
observed in the dust emission results from the inner trun- 
cation of the disk, perhaps due to planets, which does not 
effect the density distribution at larger radii. In order to 
investigate the nature of LkCa 15 disk we here adopt 
both approaches which we will be referred to as smooth 
and truncated viscous models. 

The radiation transfer equation within the disk is 
solved by adoptin g the " two-layer" approximation of |Chi- 



ang & Goldreich ( |1997[ ), in which the disk is described 



through a warm surface layer and a cooler midplane in- 
terior. Both temperatures are calculated as a function of 
the orbit al radius by iterating on t he vertical disk struc- 
ture (see Dullemond et al. 2001). The disk is in hy- 
drostatic equilibrium between the gas pressure and the 
stellar gravity, which leads to a flared geometry with the 
opening angle increasing with the distance from the cen- 
tral star. Because a flared disk geometry tends to over 
predict the far-IR disk emission in LkCa 15, we allow the 
disk scale height H p to be lower than the self-consistent 

solution H p by introducing a settling parameter \I/ < 1, 
so that H„ = ^H„. F ollowing |Mulders et alT| ( |2Q1Q| ) and 
|Espaillat etH] fl2007| ), we chose W = 0.3. 

The dust opacity is cal culated by assumin g an inter- 
stellar grain composition ( |Pollack et al.|1994 ) and a par- 
ticle size distribution between a minimum (a m i n ) and a 
maximum (a max ) size according to n(a) oc a~ q . The 
spectral slope of the LkCal5 mm-wave disk emission is 
a — 3.4, as computed from the observed flux density 
at 1.3 mm (128 mJy; see S ection 2) and 7 mm (0.44 



mJy; Rodmann et al. 
sion at these wavelengt 



. [2006]). Assuming that the emis- 
ins is optically thin and follows the 
Rayleigh- Jeans approximation, this implies a dust opac- 
ity slope /3 = 2 — a = 1.4, which is reproduced adopting 
&min = 0.5 /im, a max = 0.5 mm, and q = 3.5. For sake 
of simplicity we assume that the dust opacity is constant 
throughout the disk. 

From the derived dust density, temperature, and opac- 
ity we calculate the disk SED and synthetic disk images 
in the dust continuum at 1.3 mm. We assume that the 
center of the disk corresponds to the center of symme- 
try of the dust emission shown in panel (c) of Figure [I] 
The synthetic disk images are then Fourier transformed 
and sampled at the positions on the (u,v) plane corre- 
sponding to the CARMA observations. Best-fit values for 
the free parameters that define the disk surface density, 
and for the disk inclination and position angle are then 
derived through a x 2 minimization based o n a Markov 
Chain Monte Carlo simulation as discussed in llsella et al.l 



(2009). 



5. RESULTS 



1 | Andrews et al.| ( |2011a| adopt a different parameterization 
the similarity solution tor the surface density, which depends 
the characteristic radius r c = rt X [2(2 — 7)] 1 /( 2- ^) and on the 
normalization value E c = £ t x [2(2 - j)}^/^~ 2 \ 
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Table 1 

Best-fit parameters for the smooth and truncated models. 
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Note. — The numbers in square brackets are either fixed values (i.e., the normalization radius ro) or are the result of the model fitting. 
The inclination and position angle are expressed in degrees, the radii r*i n , r ou t, rt, and ro in AU, the disk surface densities St and Sq in g 



cm , and the disk mass in solar masses. 

5.1. Constraints on the disk surface density 

The best-fit parameters for both the smooth and the 
truncated disk models are listed in Table [I] The quoted 
uncertainties correspond to the 68.3% confidence level 
and were calculated by integrating the marginalized 
probability distribution resulting from the Monte Carl o 
fitting procedure (see Appendix A of Isella et aL||2009| ). 
In addition to the free model parameters, we list within 
the square brackets some relevant values for the disk 
structure, such as the inner and outer disk radius for 
the smooth viscous model, the outer radius and 7 for the 
viscous truncated disk model, and the normalization ra- 
dius for the power-law model. Finally, the last column 
of the table shows the total disk mass. In addition to 
the statistical errors, the surface density E t and the to- 
tal disk mass are affected by systematic errors due to 
the flux calibration and to the uncertainty on the dust 
opacity. We estimate a 10% uncertainty from flux cal- 
ibration based upon CARMA flux calibration monitor 
programs. The uncertainty on the dust opacity is more 
difficult to quantify and might vary between 5% and 20% 
once the slope of the millimeter SED is used to constrain 
the grain size dist ribution (see the discussion in Isella 
[eTatl[2009l |2010a| ). 

Figure [5] presents the comparison between the best-fit 
models and observations. Panels (a) show the real part 
of the observed correlated flux (points with error bars), 
the visibility profile of the best-fit model (solid line) and 
the residuals obtained by subtracting the model visibility 
to the observations. Panels (b) and (c) show the maps 
of the residuals obtained adopting the same weighting 
functions and color scale of panels (b) and (c) of Fig- 
ure [l] respectively. Finally, the surface density profiles 
corresponding to the best-fit models are shown in Fig- 
ure [3J along with the corresponding uncertainties. 

We find that the smooth viscous and the power-law 
models provide acceptable fits to the observations. The 
first model is characterized by 7 = —2.15, which is more 
negative, although consistent within the uncertainties, 
with the 7 value derived in 2009 by fitting only the 
smaller spatial frequencies. Both models suggest a disk 
inclination of 51° ±1° and a position angl e of 64° ±2°, 
that agrees with the values derived by Andr ews et al.| 
( |2011a[ ). By contrast, the truncated viscous model char- 
acterized by 7 = 1 does a worse job in reproducing the 
structure of the dust emission and leads to significative 
residuals both inside the continuum cavity and along the 
dusty ring (second row of Figure |2j. The comparison 
between the best-fit surface densities (Figure [3| shows 
that the first two models are characterized by rather flat 
density profiles between about 50 and 120 AU, while the 



truncated viscous model produces a much steeper de- 
crease in the surface density in this radial range. This re- 
sult therefore implies that the surface density in LkCa 15 
disk is inconsistent with the a-constant prescription at 
least in the radial range probed by our observations. 

The power-law and the smooth viscous models repro- 
duce equally well the ring observed in the dust emission. 
To understand why our observations cannot disentan- 
gle these two models, we show in Figure [4] the differ- 
ence between the correlated flux of the power-law and 
the smooth viscous model as a function of the spatial 
frequency (solid curve), along with the la and 2a noise 
level provided by our observations (blue and red boxes 
respectively). The maximum difference between the two 
models is about 6 mJy at 580 kA and is within the 2a 
noise level. The figure also shows that our ability to 
disentangle different surface density profiles is not ham- 
pered by the angular resolution (i.e., by the maximum 
spatial frequency achieved by our observations), but by 
the noise level on spatial frequencies between 100 kA and 
700 kA where the difference between the two models is 
the greatest. This is an important point to consider in 
planning futur e ob servations and will be discussed fur- 
ther in Section [6^21 

Although we cannot achieve a model-independent con- 
straint of the disk surface density, our observations pro- 
vide an upper limit of the amount of small solid particles 
present in the LkCal5 inner disk. We find that the upper 
limit is set by the best-fit smooth viscous model, since 
any additional amount of dust in the inner disk provides 
a worse fit to the data. Integrating this density profile 
between 3-42 AU, we calculate a dust mass of 7 Earth 
masses (Mq), or about 2 Mj of gas assuming a standard 
gas-to-dust ratio of 100. Additional constraints of the gas 
density come from spatially resolved observations of th e 



2007). 



12 CO and 13 CO (2-1) line emission flPietu et ah 
These observations are characterized by an angular res- 
olution between l // -2 // , and were compared to power-law 
disk models to derive the gas den sity and temperature 
radial profile. |Pietu et aL| ( |2007| find that CO might 
be present down to a radius ol 10 AU, and that a sharp 
truncation of the gaseous disk at 50 AU is excluded at 
the 7a level. These observations also suggest the pres- 
ence of about 2 Mj of gas within 42 AU from the central 
star. The amount of dust within 42 AU from the central 
star is also constrained by t he i nfrared SED modeling, 
which is discussed in Section F 



Finally, w e confirm the discrepancy found by |Pietu 
et al.| ( 2007 ) between the disk outer radius inferred from 
the dust continuum emission, i.e., about 160 AU, and 
that derived from the observations of the molecular line 



Isella A., Perez L., and Carpenter J. 



(a) 




Smooth viscous 
model 



Truncated viscous 
model 




500 1000 
Buv (klam) 



(b) 


- 1 






k + - 




/ 







(c) 







\ 


; i.. 


V. 


- to 

- \ 


+ v 


\ 

. .-'\; 
. v "\ 


m , : 


\ 

• 


. / - 

\ '' /: ~ 
\ / 




- 1 














\ en > 

\ y 






2 


1 0-1-2 




(5R.A. (arcsec) 




0.5 -0.5 

SR. A. (arcsec) 



Figure 2. Comparison between the observations and the best-fit models calculated assuming a smooth viscous (top panels), a truncated 
viscous (middle panels) and a power-law (bottom panels) surface density profile. Panels (a) show the real part of the correlated flux 
observed at 1.3 mm (square points with error bars), along with the best-fit model (solid line) and the residuals. To account for the disk 
geometry, the data have been deprojected using the values for the inclination and position angle listed in Table fl] Panels (b) and (c) show 
maps of the residuals obtained with the weighting functions used in panels (b) and (c) in Figure^ The color scale, the map size and noise 
level are also the same as in Figure^ Solid contours are plotted every 3cr, while the dotted contours in panels (c) stars at 2cr and are 
spaced by la. The dashes ellipses show the inner and outer disk radius corresponding of the best-fit models listed in Table ^ 



emission, which, depending of the observed molecule, 
varies between 550 AU ( 13 CO) and 900 AU ( 12 CO). 
Gaseous disks are found systematically larger than dusty 
disks in al most all the objects observed at high ang ular 
resolution fllsella et al.|2007| |Pietu et al.||2005[ 120071 ). In 
some cases this discrepancy can be reconciled by adopt- 
ing an exponentially tapered surface density pr ofile, as 



et al. 2008 Isella et 



the one resulting from the dis k viscous evolution ( Hughes 

aLJL 2009 |)- This model predicts that 
rthi 



dust is present up to the outer radius inferred from the 
gas but the millimeter-wave dust emission coming from 
the outer disk is too low to be detected by existing ob- 
servations. In the case of LkCa 15, we have used the 
gas emission model discussed in |Isella et al.| ( |2007[ ) to 
calculate the 12 CO J=2-l emission corresponding to the 



smooth and truncated viscous best-fit models for the 
dust continuum emission. The 12 CO surface density is 
described by Equation [I] with 7 and r t as in Table 1, 
and assuming a CO/H2 abundance of 10 -4 , as measured 
in the molecular clouds. By assuming that the 12 CO 
molecules are in thermal equilibrium with the dust, we 
find that the 12 CO emission should be confined within 
300 AU and 550 AU from the central star for the smooth 
and truncated viscous model respectively, despite the 
fact that the 12 CO density extends in practice to infinity. 
The fact that 12 CO is observed out to a radius of 900 
AU suggests therefore that the exponential taper that 
characterizes the viscous surface density is not a good 
explanation for the different radial extent of the gas and 
dust emission in LkCa 15. Other possible explanations 
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Figure 3. Surface density profiles corresponding to the smooth 
and truncated viscous best-fit models (green and blue respectively), 
and to the power-law model (red). The colored regions represent 
the 3a uncertainty interval. 
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Figure 4. Difference between the 1.3 mm correlated fluxes (solid 
curve) of the best-fit power-law and smooth viscous models. The 
blue and red triangles show the la and 2a noise level of the 
CARMA observations. The figure shows that the sensitivity of 
CARMA observations is not sufficient to distinguish between the 
two surface density parameterizations. 

are discussed in Section [631 

5.2. Sub structures and asymmetries 

Figure [2] shows the map of the residuals obtained by 
subtracting the best-fit models from the observations in 
the Fourier domain, and by adopting different weighting 
functions in the image reconstruction process. Panels (b) 
and (c) adopt the same weighting functions, use the same 
color scale, and have the same size of the maps shown in 
the panel (b) and (c) of Figure [I] 

Panel (b) clearly shows the asymmetric emission in the 
northwest side of the disk discussed in Section [31 The 
integrated flux is about 8 mJy, which corresponds to a 
dust mass between 5 and 10 Mj for a dust temperature 
between 20-40 K. 

As for the residuals along the dusty ring and inside the 
cavity, we discuss only the case of the smooth viscous 
and power-law model, since we have concluded in the 
previous section that the truncated viscous model does 
not provide as good of a fit to the observations. For 
these two models, the brightest residual in panels (c) 
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Figure 5. Radial variation of the noise in the map presented in 
panel (c) of Figure [2] The circles show the standard deviation (<x) 
of the intensity measured in concentric 0.2 /x wide annuli. The filled 
circles correspond to annuli in which the intensity distribution is 
consistent with a standard normal distribution, as expected if the 
emission is due to gaussian noise. The open circles indicate the 
radii where the intensity distribution is dominated by true emis- 
sion. The solid curve represents the theoretical noise as computed 
from the weights on the visibility data. 



is identified with the letter A and has peak flux of 2.4 
mJy. It is located along the major axis of the disk at 
the physical distance of 140 AU from the center, and 
therefore very close, or just outside, the outer disk radius 
inferred from our disk modeling, which is represented by 
the large dashed ellipses in Figure [2] At this distance 
from the central star, the dust is expected to have a 
temperature of about 40 K, which leads to a mass of 
dust of about 1.3 Mj. 

The significance of this residual depends critically on 
knowing the value of the noise in our map and how it 
varies with the distance from the center due to the an- 
tenna primary beam attenuation. A theoretical value for 
the noise can be obtained from the weights on each visi- 
bility point, which are measured at the time of the obser- 
vations and depends on the effective are a of the telescope 



and t he system temperature (see, e.g, Thompson et al 



1991). However, if the observed source is bright, this es 



timate does not take into account possible effects due to 
the image deconvolution, such as the presence of residual 
side lobes in the final map, and dynamic range limita- 
tions due to the atmospheric decorrelation (i.e., seeing) 
and calibration errors. 

We estimated the noise in the images as follow. We 
divide the intensity map in concentric annuli character- 
ized by a width of 0.2". For each annulus we calculate 
the mean intensity and the standard deviation a. If the 
emission is due only to gaussian noise, then the mean is 
and the standard deviation of the intensity is equal to 
the noise. For each annulus, we check that the distribu- 
tion is gaussian by creating a histogram for the intensity 
and fitting a normal distribution. Otherwise, if the an- 
nulus contains true emission, the intensity distribution 
will deviate from the normal profile. In this case a will 
give an upper estimate of the noise. 

Figure [5] shows the radial variation of the standard de- 
viation ofthe intensity, a. At radii larger than 6" (filled 
circles), the intensity distribution is well described by a 
standard normal function; a increases with the radius as 
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the inverse of the antenna primary beam profile and it 
is equal to the theoretical noise indicated by the solid 
curve. Within 6", the intensity is a combination of gaus- 
sian noise and true emission and its distribution deviates 
from the standard normal profile. In this region, the 
noise in the map may assume values between the theo- 
retical noise and a. This implies a noise upper limit of 
0.47 mJy/beam at the position of the residual A (i.e., 
1"), which leads to a signal-to- noise ratio of 5.1, given 
the measured peak flux of 2.4 mJy. 

The significance level of residual A can be estimated 
by multiplying the normal probability corresponding to 
a 5.1<r residual (3.4 x 10 _T ) by the number of inde- 
pendent pixels in our map. If each resolution element 
(0.21" x 0.19") is accounted as independent pixel, then 
the number of independent pixels within the 1.3 mm 
CARMA field of view (27") is about 18000, and the prob- 
ability that the residual might be due to random noise is 
smaller than 0.6%, or about 2.8<r. This result suggests 
that the residual might trace a real compact over-density 
in the dust distribution with respect to the best-fit sur- 
face density profile, but more sensitive observations are 
required to investigate the nature of this emission. 

All the other residuals in panel (c) of Figure [2] have flux 
within ±4.5 times the noise level. The probability that 
any one of these residuals could originate from random 
noise is greater than 10%. 

5.3. Understanding the SED 

Figure [6] compares the observed SED for LkCa 15 with 
the SED inferred from the best-fit models to the millime- 
ter data. Both the smooth viscous (upper panel) and the 
power-law (lower panel) provide a good agreement with 
the observations at millimeter wavelengths. At infrared 
wavelengths, the model SED strongly depends on the 
disk surface density profile. In the case of the smooth 
viscous model, the disk is optically thick to the stellar 
radiation starting from 3 AU, and the agreement with 
the observation extends to about 20 jam. At shorter 
wavelengths the model SED underestimates the observed 
emission from the disk. This is due to the fact that the 
surface density within a few AU from the central star 
is lower than 10 -6 g cm -2 and the near-infrared disk 
emission is negligible with respect to the stellar photo- 
sphere. In the case of the power-law model, the predicted 
SED agrees with the observation only at wavelengths 
longer than 50 jam due to the cut-off of the dust den- 
sity at 42.5 AU. Independently on the adopted model, 
the near-infrared emission observed toward LkCa 15 re- 
quires the presence of material with temperature between 
1000-2000 K in excess of what predicted by our models. 

For a power-law surface density profile truncated at 
42.5 AU, the near infrared emission ca n be explained 
by tw o different models, as discussed in Mul ders et al. 
(2010). In the first case, the emission comes from a small 



optically thick disk that extends from the dust s ublima- 
tion radius C^O-l AU ) up to about 1 AU (see also Espail- 
lat et al.|2007 |2008 ). For typical infrared dust opacities, 
this model requires a minimum dust mass of 3 x 10 -5 
M®. Alternatively, the LkCal5's infrared SED can be 
fitted with an optically thin region that extends from 0.1 
AU to 5 AU. In this case, the mass of dust would be 
3 x 10 -6 M®. The optically thin region and the optically 
thick disk would produce a 1.3 mm flux below 10 -6 mJy 
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Figure 6. SED of LkCa 15. Filled squares present phot omet- 
ric data from the literature (Kenyon & Hartmann 1995 Hartmann 
et al. 2005 Espaillat et al. 2007 Kebull et al. 2U1U Kitamura et al 
2002| jDutrey et al. | l996| | Andrews fe Williams|2005| |2007| |Oberg 
|et al.|20l0[ and from our new CAKMA observations. Infrared pho- 
tometry at 9 jLim, 19.6 fim, and 90 /xm is from the Akari archive, 
and the infrared spectrum from 5 /im to 14 /im is from archival 
Spitzer IRS observations. The solid curves show the SED of the 
best-fit models for the smooth (top panel) and the power-law (bot- 
tom panel) surface density parameterization. The dashed curves 
indicate the SED required to reproduce the observed near- and 
mid-infrared excess. As discussed in the text, this additional com- 
ponent can arise either from a narrow ring of optically thick dust, 
or from a more extended optically thin region. 

and would therefore have a negligible effect on the overall 
disk millimeter emission. The same two models for the 
near-infrared emission can be applied to the case of the 
smooth viscous surface density profile. However in this 
case the outer disk is brighter in the mid-infrared than 
in the power-law case, and therefore the region emitting 
the infrared emission has to be smaller. 

6. DISCUSSION 
6.1. On the origin of the continuum cavity 

The results presented so far suggest at least two pos- 
sible scenarios for the radial distribution of the circum- 
stellar material around LkCa 15. 

In the first case, most of the circumstellar dust is con- 
fined in a sharply truncated ring extending between 42- 
120 AU. Since no discontinuity is observed in the gas 
density at the edges of this dusty ring, this scenario re- 
quires a drop of the dust-to-gas ratio, or of the dust opac- 
ity, outside the ring. As discussed in Section [I] a sharp 
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internal disk truncation can be the result of dynamical 
interactions with o ne, or more, giant planets (see, e.g., 
Bryden et al.||1999 ) or be due to the disk photoev apora- 



tion caused by the stellar radiation field (see, e.g., Gorti 
|et al.|[2QQ9| ). This latter model suffer however two se- 
vere problems. First, it predicts a truncation in both 
the gase ous and dusty disk, which, as discussed in Sec- 



tion [5J-] it is not consistent with spatially resolved obser- 
vations of the mm-wave CO emission. Second, it predicts 
that the accretion of material on the central star stops 
as soon as the cavity is opened, while the mass accretion 
rate meas ured toward LkCa 15 is about M = 3 x 10 -9 
M yr- 1 dHartmann et al . |1998| > . 

Whether large (> 20 AU) dust depleted cavities can be 
produced by the dynamica l interaction with giant planets 
is still subject of debate. Zhu et al. (2011) showed that 
for large values of the viscosity parameter (a > 0.01), a 
system with 3 or 4 Jupiter mass planets orbiting between 
5 and 20 AU can lead to a largely depleted gap in the 
disk surface density. However, the mass accretion rate on 
the central star is predicted to drop below 5 x 10 -10 M 
yr -1 in less that 1 Myr after the formation of the planets. 
Furthermore, a discontinuity in the gaseous disk surface 
density will not be consistent with the CO observations. 
Alternatively, a lower viscosity (a < 5 x 10 -3 ) would lead 
to larger values of M and less depleted cavities. While 
this latter model might be consistent with the observa- 
tion of the CO emission, it fails to explain the lack of 
mid-infrared dust continuum emission observed toward 
LkCa 15. To reconcile the measured values of M with 
the i nfrared SEP in th e framework of planet- disk interac- 
tion, |Zhu et ah] (2011 ) argue that low values of a have to 
be coupled with a decrease by at least one order of mag- 
nitude of the near- and mid-infrared dust optical depth 
in the region partially depleted by the planets. A drop 
of this magnitude in the dust opacity can be achieved 
through the coagulation of small, micron sized grain into 
larger bodies. In the particular case of LkCa 15, we find 
that the dust opacity at 10 /im will decrease from about 
3 cm 2 g _1 to 0.2 cm 2 g _1 by increasing the maximum 
grain size from the adopted value of 0.5 mm (Section [4| 
to 10 cm. It is therefore possible that the actual structure 
of the LkCa 15 inner disk is the result of an earlier phase 
characterized by the formation of large peb bles and boul- 
ders i n the innermost disk regions (see, e.g., Brauer et al. 
2008 ) , followed by the formation of giant planets. 11 that 
happened in a low viscosity environment, the gas density 
within 40 AU would be high enough to explain both the 
measured mass accretion rate and the CO observations, 
and perhaps to allow the planetary cores to accrete the 
gas required to reach the Jupiter mass. 

This scenario is compatible with the recent discovery 
of a possible proto-planet orbiting at about 16 AU from 
the central star discovered through aperture masking ob- 
servations in the K' [~2.l fim) and V (~ 3.8/im) bands 



flKraus fc Ireland|2011[ ). The K' magnitude and K' - V 
color ol the planet are consistent with a temperature of 
1500 K and a mass of 6 Mj, or smaller. In addition the V 
observations trace extended emission interpreted as ev- 
idence of circumplanetary material. However, following 
theoretical models for the planet disk interactions, this 
planet would be able to open only a relatively small gap 
in the dust distribution between 12-20 AU. To explain 
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Figure 7. Dependence of the slope of the dust opacity (3 mea- 
sured between 1.3 and 3 mm (solid curve) and of the dust opacity k 
at 1.3 mm (dashed curve) on the maximum grain size, for a grain 
size distribution n(a) oc a~ q , where q = 3.5 and the minimum 
grain size is 5 x 10 — 6 mm. 

the cavity radius of 42 AU, at least another planet clear- 
ing the disk within 12 AU, and other two planets clearing 
the region between 20 and 42 AU wou ld be required (see, 
e.g., |Dodson- Robinson fc S alyk 2011). 

In an alternative scenario, which does not require mul- 
tiple planets, the disk surface density is described by the 
smooth viscous model shown in Figure [3j In this case, 
the gas density decreases within 80 AU but it remains 
high enough to explain the small inner disk radius de- 
rived from the CO observations. This situation is similar 
to that proposed for MWC 758, where the 12 CO emission 
is centrally peaked while the dust emission traces a par- 
tially ^dus^_de_2leted_d^k within 70 AU from the central 
star ( |Isella et al.|2010b| ) . It is important to note that the 
smooth best-tit model was derived in the assumption of 
constant dust opacity throughout the disk. By relaxing 
this assumption, we can interpret our results in terms of 
radial variation of the dust opacity, due, for example, to 
grain growth. Radial variations in the grain size distri- 
bution are indeed predicted by theo retical models for the 
dust evolution ( |Brauer et al.||2008| ) and can be observed 
by measuring radial variations ol the mm-wave opacity 
slope /3 through spatially resolv ed multi- wavelengths ob- 
servations of the dust emission ( |Isella et al.| p 
In the case of LkCa 1 



loteau et al. 



|2010a 



Guil- 



com- 



2011]). In the case ol LkCa 15, the 
parison ol O and 3 mm observations suggests that j3 
incr eases from 0.7-1.0 to 1.5 -1.7 between 30 AU and 150 
AU flGuilloteau et al.pQll]) . As shown in Figure [F| this 
variation would correspond to a decrease of the maxi- 
mum grain size from about 8-30 mm to 0.3-3 mm, or to 
an increase of about a factor of 2 in the dust opacity at 
1.3 mm. The observed variation in the grain size might 
therefore explain the observed cavity in the 1.3 mm dust 
emission without invoking multiple giant planets, but 
more sensitive observations at longer wavelengths are re- 
quired to investigate the presence of large grains in the 
1.3 mm continuum cavity. 

6.2. Investigating the origin of transition disks with 
ALMA 

Several experiments can be set up to investigate the 
origin of the inner cavities in transition disks by taking 



10 



Isella A., Perez L., and Carpenter J. 



Smooth viscous model 



o o 



CD 
Q 



— I — I — I — I — | — I — I — I — I — | — h 

Power— law model 



ALMA Simulation 




-| — I — I — I — I — | — I — I — I — I — | — h 

ALMA Simulation 






6R.A. (arcsec) 





6R.A. (arcsec) 



Figure 8. Simulated ALMA observations of 0.44 mm continuum 
emission corresponding to the smooth viscous (top panels) and 
power-law (bottom panels) best-fit models discussed in Section [5?T] 
The panels on the left show the model while the simulated observa- 
tions are shown on the right column. The simulations assume a 30 
min on source integration and the default weather conditions for 
Band 9. ALMA will determine whether the inner edge of the disk 
is sharply truncated (possibly due to dynamical interactions with 
a planet), or varies smoothly with radius (possibly due to grain 
growth in the inner disk). 

advantage of the superb imaging capabilities of the Ata- 
cama Large Millimeter Array (ALMA). Here we explore 
the possibility to use relatively low angular resolution 
observations, achievable during ALMA early science, to 
constrain the profile of the dust distribution at the edge 
of the dust continuum cavity. These observations can 
be used to disentangle the two possible scenarios dis- 
cussed abov e. As shown in Figure [4j and discussed in 
Section [57T| CARMA observations fail to distinguish be- 
tween a sharp and smooth edge of the continuum cavity 
because they do not provide enough sensitivity on the 
spatial frequencies where the difference between the two 
models is the largest, namely between 100 and 700 kA. 
On the other hand, Figure [8] shows that constraining the 
radial density profile at the edge of the continuum cav- 
ity is within the capabilities of ALMA in the early sci- 
ence phase. The left-most panels of Figure [5] show the 
model of the dust emission corresponding to the best-fit 
solutions for the smooth viscous a nd p ower- law surface 
density profile discussed in Section |5.1| The maps in the 
right-most panels represent ALMA band 9 (A = 0.44/im) 
simulated observations obtained using the array specifi- 
cations for Cycle observations and 30 min on source 
integration. The structural differences between the two 
models are evident in these images: the smoothly vary- 
ing E(r) shows trace amounts of emission in the inner 
disk that is detectable with ALMA, while the inner re- 
gion in the truncated disk model is devoid of emission 
and resolved. More quantitatively, Figure [9] shows the 
difference between the visibility profiles of the two differ- 
ent models as a function of the spatial frequency, where 
the colored region indicate the 3a measurement uncer- 
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Figure 9. Difference in the visibility profile between the power- 
law and smooth viscous best-fit models (solid curve), compared to 
the 3cr sensitivity of ALMA early science observations at 0.44 mm 
(shaded region). This figure demonstrates that ALMA observa- 
tions can reliably distinguish between two plausible geometries for 
transitions disks to identify the mechanism likely responsible for 
inner cavities. 

tainties from the simulated observations. At the wave- 
length adopted for the simulation, the two models can be 
disentangled at high significance level on almost all spa- 
tial frequencies shorter than 600 kA, which correspond to 
baselines shorter than 270 m. 

6.3. The LkCa 15 outer disk structure 

The discrepancy in the outer radius of the dusty and 
gaseous dis k re mains an unsolved issue since, as discussed 
in Section |5.1[ it cannot be reconciled by adopting an 
exponentially tapered surface density profile. A caveat 
to this analysis is that our observations constrain the 
dust density only between 40-120 AU. The extrapolation 
of E to larger radii relies therefore on the assumption 
that the prescription of Equation [I] is valid across the all 
disk radial extent. 

To test this assumption we compare our results with 
the gas densit y inferred from the analysis of the 13 CO (1- 
0) emission by |Pietu et aT| ( |2QQ7| ). In LkCa 15, the 13 CO 
(1-0) line is nearly optically thin and traces the density 
in the disk mid-plane. By comparing the ob servations 
to a power-law disk model, [Pietu et aL (2007) find that 
the 13 CO surface density decreases as r -3 / 2 from about 
1.4 x 10 16 cm" 2 at 100 AU to 1.2 x 10 15 cm" 2 at 500 AU. 
If the 13 CO abundance is constant throughout the disk, 
then the decrease of the gas density is less steep than 
what is predicted by an exponentially tapered profile. 

To check whether the disk surface density profile de- 
rived from CO observations is consistent with the dust 
continuum emission, we show in Figure [To] simulated 
CARMA observations of the 1.3 mm dust emission cor- 
responding to a dust surface density E^(r) = 0.13 x 
(r/100AU)" 3/2 that extends between 42.5 AU and 550 
AU. The surface density normalization is derived in or- 
der to reproduce the peak flux measured in our map. To 
facilitate the comparison with the observations, the sim- 
ulated emission is observed at the same resolution and is 
affected by the same noise as in panel (c) of Figure fl] 

We find that only the dust emission within about zOO 
AU would have been clearly detected while the flux aris- 
ing from larger radii would have been lost in the noise. 
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Figure 10. Simulated observations of the 1.3 mm dust continuum 
emission corresponding to the surface density profile derived from 
the a nalysis of the optically thin 13 CO line emission ( |Pietu et al.| 
|2007| >. The simulation have been obtained by assuming the same 
uv coverage and noise level as in the CARMA observations. The 
image as been made by adopting uniform weighting to achieve a 
resolution of 0.21" xO.19", as in panel (c) of Figure^ Contours 
are plotted every 3cr, where the la noise level is 0.43 mJy/beam. 

The limited sensitivity of the continuum observations 
seems therefore to be the major responsible for the dif- 
ferent radial extent of the dust and gas emission. Nev- 
ertheless, the fact that the simulated emission is about 
60% more extended than what is observed suggests that 
either the dust-to-gas ratio or the dust opacity decreases 
at about 130 AU from the central star. 

In Section |6.1| we suggested that grain growth is re- 
quired to explain the structure of the innermost disk re- 
gions independently from the presence of planets. Here 
we argue that grain growth might also provide a natural 
explanation for the smaller extent of the dust contin- 
uum emission. Grain growth models indeed predict that 
the outermost disk regions should be populated by sub- 
micron sized grains, which have a 1.3 mm dust opacity 
almost a factor 10 smaller than the mean dust opacity 
adopted in our disk modeling (Figure [7]). In the assump- 
tion that particles grow through mutual collisions and 
stick together by van der Walls forces, the maximum 
grain size a m , ax of the grain size distribution can be ex- 
pressed as ( |Birnstiel et al.|2010| ) 



■Xr) 



2S(r) 



nap s c s (r) 2 



(2) 



where a is the viscosity parameter, Uf is the dust frag- 
mentation velocity, E is the disk surface density, c s is the 
sound speed , and p s the dust density. Figure [TT] shows 
the radial profile of a max , and of the corresponding dust 
opacity at 1.3 mm, calculated by assuming p s = 2 g cm 2 , 
a = 10 -3 , Uf = 1 m s _1 . The surface density is that de- 
rived from the CO observations as discuss above, and the 
sound speed c s is self-consistently calculated by assuming 
hydrostatic equilibrium (see, Appendix E in |Isella et al.| 
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Figure 11. Top panel: Dependence of the maximum grain size 
CLmax (solid curve) and of the dust opacity k at 1.3 mm (dashed 
curve) on the orbital radius .calculated assuming the grain growth 
model discussed in Section [6l Bottom panel: Simulated observa- 
tions of the 1.3 mm dust continuum emission corresponding to the 
surface density profil e derived from the a nalysis of the optically thin 
13 CO line emission ( Pietu et al.||2007| ), accounting for the radial 
variation of the dust opacity due to grain growth. The decrease in 
the dust opacity between 100-200 AU results in a smaller effective 
radius for the dust continuum emission, compared to the ca se in 
which the dust opacity is constant throughout the disk (FigurefTo]). 



2009). Although the exact values of a max depend on the 
assumptions on a and w/, we find in general that a ra- 
dial decrease of the maximum grain size corresponds to 
a drop in the 1.3 mm dust opacity. In particular, the de- 
crease in a max from 0.4 mm to 0.1 mm between 100 and 
300 AU corresponds to a factor 5 decrease in the dust 
opacity. Given the sensitivity of our observations, the 
drop in the opacity would in practice lead to a smaller 
radius of the continuum emission, making the simulation 
much alike the observations (bottom panel of Figure 11 ). 

This result suggests that the small disk radius mea- 
sured in the dust continuum emission might be the re- 
sult of a combined decrease in dust opacity due to radial 
variations in the grain size distribution and the finite sen- 
sitivity of our observations. Investigating this hypothesis 
requires high-angular resolution multi-wavelength obser- 
vations of the dust continuum emission. 

7. SUMMARY 

We report new CARMA 1.3 mm dust continuum obser- 
vations toward LkCa 15 which resolve the disk on spatial 
scales of 30 AU. We find that 90% of the observed 1.3 
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mm dust emission arises from an azimuthally symmetric 
ring that contains about 5 x 10 -4 M of dust. The re- 
maining 10% of the emission comes from a low-surface 
brightness tail that extends toward the northwest out to 
a radius of about 300 AU. 

To investigate the origin of the inner cavity observed 
in the dust emission we constrain the radial profile of 
the dust density at the edge of the cavity and across the 
dusty ring. The comparison with theoretical disk models 
shows that the surface density is rather flat between 40 
and 120 AU, and that the drop of dust emission is consis- 
tent with either a sharp truncation of the dusty disk at 
42.5 AU, or with a smooth inward decreases of the dust 
density between 3-85 AU. Within 40 AU, the observa- 
tions constrain the amount of dust between 10 -6 and 7 
, where the minimum and maximum limits are set by 
the near infrared SED modeling and by the millimeter- 
wave observations of the d ust emission respectively . 

The recent discovery by Kraus & Ireland (2011) of a 
possible giant planet orbiting at 16 AU around LkCa 15 
suggests that the dust distribution might be shaped by 
the dynamical clearing of this planet. However, at least 
3 other Jupiter size planets in addition to the one discov- 
ered by Kraus & Ireland would be required t o explain a 



sharp truncation in the dusty disk at 42 AU ( Zhu et aL 

Salyk|[201l| . In addition, the' 
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presence ol planets should be coupled with low disk vis- 
cosity (a < 5 x 10 -3 ) and with a significant depletion of 
micron-size grains, to explain both the presence of gas 
within the dust deplete d cavity and the l ack of near- and 
mid-infrared emission (Zhu et al. 2011). Alternatively, 
our results are consistent with a smooth inward decrease 
in the dust opacity caused by the coagulation of small 
grains into larger bodies, as predicted by grain growth 
models. In this case, the dynamical clearing by the dis- 
covered planet would play a minor role in shaping the 
disk emission and no additional planets are necessary. 
We show that early science ALMA observations will be 
able to disentangle these two scenarios by constraining 
the slope of the dust density at the edge of the continuum 
cavity. 

Finally, we find that 90% of the dust emission is con- 
fined within a radius of about 120 AU, while the gas 
emission is observed up to a radius of 900 AU. We find 



that this discrepancy cannot be reconciled using a expo- 
nentially tapered surface density profile but suggests a 
drop in the dust-to-gas ratio, or in the dust opacity, at 
radii larger than 120 AU. Using the grain growth model 
from Birnstiel et al. ( 2010[ ), we show that a decline in 
the maximum gram size with radius will reduce the dust 
opacity by a factor of 5 between 100 and 300 AU, with a 
corresponding decrease in the dust continuum emission. 

Note added in proof: 

Since the work described in this paper was completed 
and submitted for publication, 850 [im observations 
of the dust continuum emission toward LkCa 15 that 
ac hieve a resolution of abou t 0.25" have been published 
by | Andrews et al.| ( |2011b|) . These observations show 
that the continuum emission arises from a ring with 
the same properties of that discussed in this paper. In 
particular, they also find that the ring requires a rather 
fiat surface density profile between about 40 and 120 AU. 
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APPENDIX 
STELLAR PROPER MOTION 



Th e proper motion for LkCa 15 has bee n measured by Roeeser et al. J2010[) comparing USN O-Bl-0 (Monet et al. 
2003) and Two Micron All Sky Survey data ( Skrutskie et ar|2t)06| ) , and by Zachaxias et al. (|2009| US Naval O bservatory 
Twin Astrograph project). The first measurement leads to [i a cos S = 9.1 ± 2.1 mas yr and [i§ = —13.8 ± 2.1 mas 
yr _1 . The second to [i a cos5 = 9.6 ± 2.2 mas yr -1 and [is = —13.3 ± 2.2. The J2000 c oordinates for epoch 2000 



for LkCa 15 corrected for the proper mot ion are RA=69.824139 ° and Dec=+22. 350945° ( |Roeeser et aL||2010| ), and 
RA=69.8241427° and Dec=+22. 3509412° (ISkrutskie et al.||2006|). The difference between the two catalogs is about 



13 mas, and is cons istent with the quoted astrometric errors (i.e, 15 mas in both right ascension and declination in 
Roeeser et al.|201Q , and 30 mas in Zacharias et al.|[2009 respectively). Our CARMA observation have been corrected 
by assuming a mean proper motion of [i a cosS = 9.35 mas yr -1 and [is — —13.55 mas yr -1 . The J2000 coordinates 
corrected for the proper motion at the epoch of the observations are listed in Table [2] 
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